Spin Rotation Invariant Spin Triplet Superconducting liquids 
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Spin ordering and its effect on the low energy quasiparticles in a p-wave superconducting fluid 
are investigated. We study the properties of a new 2D quantum spin triplet superconducting liq- 
uid where the ground state is spin rotation invariant. In quantum spin disordered cases, zero 
energy skyrmions are topologically stable because of a fermionic zero mode in space-time monopoles 
which serve as a quantum protectorate. The low energy quasiparticles are bound states of the bare 
Bogolubov- De Gennes (BdeG) quasiparticles and skyrmions, which are charge neutral bosons at 
the zero energy limit. Further more, the spin collective excitations are fractionalized ones carrying 
a half spin and obeying fermionic statistics. In thermally spin disordered superconducting states, a 
bare BdeG has an infinity energy and each two BdeG quasiparticles form spin zero bound states. 



I. INTRODUCTION 

Quantum number fractionalization has been a fasci- 
nating subject. The curiosity about this subject has been 
widely shared by mathematical physicists, field theorists 
and condensed matter physicists |l|-|l(i|] . The earliest ex- 
ample of the quantum number fractionalization is per- 
haps the angular momentum fractionalization of an elec- 
tron when confined with a Dirac monopole Q. As now 
known to many of us, the electron carries a half inte- 
ger angular momentum in this case, by contrast to what 
we've learned in atomic physics. 

Zero modes in certain topological excitations are vi- 
tal to quantum number fractionalization. This was first 
appreciated in a ID field theory model where it was ob- 
served that fermionic excitations in a charge conjugate 
model can carry one half of the elementary quantum 
number H. A powerful method developed later on pro- 
vides further insights into the nature of the problem || . 
The authors of Jq] showed that in a continuum theory, 
the quantum number of a topological excitation is a con- 
tinuous function of the parameters in the theory and can 
be an arbitrary fraction in the absence of a charge conju- 
gation symmetry j^] . These ideas were also deeply rooted 
in the beautiful theories on the index of the Dirac oper- 
ators which were put forward both before and after the 
discoveries in field theories [I4I . The mechanism leading 
to the fractionalization in these cases is the spectral flow 
in the presence of a singular gauge transformation. The 
spectral flow, especially the level crossing in the field the- 
ory furthermore results in remarkable ideas of Wittcn's 
global anomaly, the Skyrmeon model of nucleaon etc [O . 

The statistics of these exotic objects in spatial dimen- 
sions higher than one was investigated at the same pe- 
riod of time. A spinlcss bosonic particle bound with a 
monopole has half spin and behaves as a fermion and vice 
verus In this situation, statistics transmutation 

takes place which further complicates the issue. The au- 
thors of [ fig] first introduced 9- statistics in two-dimension 
in their tour de force analysis of two-body wavefunctions 
with arbitrary boundary conditions. Later on, the pos- 



sibility of having anyons in 2D was discovered, which 
can be achieved by a flux attachment and which now be- 
comes a very powerful icon fl6|| . In an interesting paper, 
Wilczek and Zee further emphasised that the spin and 
statistics of a spin texture can be related to the link- 
ing number of two trajectories in the target space; the 
linking number is identical to a Hopf invariant of a rotat- 
ing skyrmion and results in anyonic skyrmion excitations 

The interest on the quantum number fractionalization 
in condensed matter systems was very much prompted 
by the remarkable discoveries of charge one but spinless 
quasiparticles and spin one but chargeless quasiparticles 
in one dimensional polyacetylene Q . The experiment dis- 
coveries of fractional quantum Hall effects in 2DEG and 
the theory on FQH incompressible liquids with quasipar- 
ticles carrying one third of an electron charge fur- 
ther stimulated investigations on this subject. This phe- 
nomenon inpired many theoretical works on fractional 
statistics. And the later theoretical investigations on 
the mechanism of high temperature superconductivities, 
where the important concept of spin-charge separation 
was introduced were once more focusd on this is- 

sue. In the backdrop of FQHEs in 2DEG and high tem- 
perature superconductivities in cuprates, two intensively 
studied subjects in condensed matter physics in past two 
decades, many efforts have been made to achieve a better 
understanding on the issue of quantum number fraction- 
alization gf-§§]. 

The theory on one dimension polyacetylene was es- 
tablished more than two decades ago || and is worth a 
detailed discussion here because the quasiparticles in 2d 
p-wave superconductors here reminisce, in terms of spir- 
its, the mid-gap states discovered in that paper. There 
are at least two essential ingredients in those one dimen- 
sional systems which make the fractionalized excitations 
possible. And we want to emphasis both of them. 

In one dimensional polyacetylene of CB-polymer, each 
carbon atom has to share with the neighboring atoms 
three a electrons and one 7r electron. The a electrons 
form covalent bonds and only ir electrons are conducting 
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and responsible for the electronic properties of the poly- 
acetylene. So the conduction band is half filled which 
is subjected to a Peierls instability. A lattice distortion 
takes place; carbon atoms at the 2n + 1-th sites along 
the chain are displaced by a 0.2A in one direction while 
atoms at the 2n-th sites are displaced towards an oppo- 
site direction. A charge density wave with a period of 
two lattice constants, 2a is formed. Simultaneously, a 
Peierls gap opens up right at the fermi surface leaving 
behind an incompressible electron liquid where the va- 
lence band is completely filled and separated from the 
empty conduction band by the Peierls gap. 

The ground state obviously has twofold degeneracy 
when the Peierls instability takes place. In a dimerized 
picture, two degenerate configurations A and B can be 
transformed into each other by a displacement of a lattice 
constant. This suggests kink- like solitons in the poly- 
acetylene, which are the domain-walls separating A and 
B configurations. A careful calculation shows that the 
size of these kinks can be rather small, of a few lattice 
constants implying a small mass for the topological ex- 
citations. The analysis carried out by Su, Heeger and 
Schrieffer convincingly demonstrates that when a kink 
like soliton is inserted into the chain, the fermionic spec- 
trum which is otherwise fully gapped, starts to develop 
a mid-gap state. Electrons of each spin specie in the va- 
lence band contribute a half mid-gap state so that there is 
precisely single electronic state with two fold spin degen- 
eracy appearing at the fermion surface. The existence 
of topological excitations, particularly that of the zero 
modes hosted in these excitations is the first ingredient 
we wish to emphasis. 

When the mid-gap states are doubly occupied or 
empty, the corresponding quasiparticles are spinless but 
carry a unit negative or positive charge. When the mid- 
gap states are singly occupied, the quasiparticles carry 
no charge but with spin one half. The charge and spin 
conservation are guaranteed by the conservation of topo- 
logical winding numbers. For a give finite system, one 
always has to insert a kink-anti kink pair. When the 
kink carries spin one half (up) but zero charge, the anti- 
kink has to carry spin down and zero charge. However, 
the kink can also carry a unit positive charge but zero 
spin while the antikink carry a unit negative charge and 
zero spin. 

It is also obvious that quasiparticles are weakly in- 
teracting and well-defined only when the topological ob- 
jects, in this case, the kinks and anti-kinks don't have 
long range interactions. This is in fact a necessary con- 
dition for fractionalization to take place and is the sec- 
ond ingredient we are turning to. In a mean field ap- 
proximation which turns out to be a remarkably good 
approximation in this case, the two fold degeneracy in- 
herited in the Peierls instability ensures the absence of a 
long range interaction and a liberation of kink-anti kink 
pairs. However, this is not so when the degeneracy is 



absent. For instance in polythiophene, a cousin of poly- 
acetylene, the two-fold degeneracy is lifted. The kink 
and anti-kink interact with a potential which grows lin- 
early as a function of distance between them. This is 
a sign of confinement of kinks as well as the fractional- 
ized quasiparticles discussed above. In a most general 
situation, the interaction between fractionalized objects 
can be highly nontrivial. Especially in many magnetic 
systems, including the one we are going to address in 
this paper, fractionalization takes place only when a spe- 
cial mechanism of quantum protection exists. And this 
quantum protectorate functions only under certain topo- 
logical enviornments, mostly through conserving certain 
topological charges or ordering of topological fields. 

To demonstrate the idea of fractionalization in mag- 
netic systems, we turn to a brief discussion of ID antifer- 
romagnetic spin chain, which has been a subject of ex- 
tensive investigations and is relatively better understood 
than other higher dimension systems. For this purpose, 
we address this issue in the context of the AKLT formu- 
lation of the problem so that a close analogy can be drew 
between the spin chain and one dimension charge con- 
jugate polymer discussed above [^9). For spin one half 
chain with Majumdar- Ghosh Hamiltonian 

H = J2 P 3/2(Si + S i+1 + S i+2 ) (1) 

i 

where P3/2 is the projection operator for a spin-3/2 state, 
the ground state can be viewed as a collection of spin 
singlet pairs formed between adjacent sites. If each site 
emits one monomer standing for a one-half spin, the 
ground state can be represented by a periodical array 
of dimcrs connecting each pair of two adjacent sites, or a 
valence-bond crystal so to speak, with a period of two lat- 
tice constants. The ground state has twofold degeneracy, 
the configurations of which (A and B) can be obtained 
by displacing one array of dimers by a lattice constant. 
The excitations are fully gapped and the spin system is 
incompressible. 

A spin one half excitation can be created by break- 
ing a dimer into two monomers and moving one of the 
monomers to the infinity while keeping the other one 
fixed in the bulk of the spin chain. So again, a monomer 
separating the A and B configurations corresponds to a 
kink; this time, one half spin is hosted in the kink. What 
one can easily show is that kinks are liberated in this case 
and spin- 1/2 are free excitations because of the two fold 
degeneracy in the problem. This is not surprising at all 
given that we started with spin-1/2 chain. 

Let us now turn to a less trivial example examined by 
AKLT. Consider a spin-1 chain with the following Hamil- 
tonian 

H = 2^P 2 (S, + S l+1 ) -(/?+i)(S,-S m ) 2 . (2) 

i 
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At (3 — —1/3, the ground state again is a valence bond 
crystal; two monomers are emitted from each site and 
terminated at two different neighboring sites. Each two 
adjacent sites are connected by one dimer forming a VBC 
with a period of one lattice constant. The ground state 
is nondegenerate and the excitations are fully gapped. 

Creation of spin 1/2 excitations, which carry only half 
of the spin of the underlying particles at each site, in- 
volves breaking a dimer and sending one of the monomer 
to infinity as emphasised above. But, because the ground 
state is nondegenerate, this procedure unavoidly provoks 
some higher energy configurations, say a string of dimer- 
ized pairs(see Fig.l). In a dimerized pair configuration, 
two monomers emitted from each site always end at a 
same site so that a singlet is formed between two ad- 
jacent sites. A string of dimerized pairs emitted from 
a spin 1/2 excitation, or a monomer can be terminated 
either at the infinity or at another spin 1/2 excitation. 
This suggests a confinement of spin- 1/2 excitations in 
spin chains with underlying spin equal to one and we are 
not able to create a free excitation which carries half of 
the underlying spin. This is again not surprising. 
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Fig.l Creation of spin-1/2 excitations in ID VBC states 
in a spin one half chain a) and a spin one chain b). In b), 
to create two spin 1/2 excitations (monomers), a dimer- 
ized string (down figure) has to be introduced to the 
ground state (up figure). 

What is important and interesting is that the energy 
difference of a dimerized pair configuration and the VBC 
can be adjusted by varying (3. As (3 approaches a large 
positive value, the energies of dimerized configurations 
are lowered and the dimerized states are stablized. A 
large negative value of j3 on the other hand stablizes the 
VBC state. It is therefore expected that these two con- 
figurations have to be degenerated at a certain point, 
mostly likely at /3 = 1. When this happens, the confin- 
ing potential disappears encouraging a spin fractionalizcd 
excitation with half spin according our argument. This 
oversimplied argument, which doesn't take into account 
fluctuations etc nevertheless demonstrates the key idea 
behind fractionalization. And a surprise doesn't come 
at no cost. A spin-1/2 excitation in a spin chain of un- 
derlying spin equal to 1 can be observed only at certain 
points {(3 = 1 for example) out of the whole parameter 



space. This time we are also lucky enough to get the 
right answer as far as the spin of the excitations is con- 
cerned regardless the crude approximation in our reason- 
ing. More sophisticated calculation using Bethe ansatz 
solutions does confirm the existence of these objects [ |30| . 
This is further supported by k — 2 conformal field theory 
calculations pjj . 

The picture just demonstrated in this short discussion 
plays a central role in the phenomena of fractionaliza- 
tion. It is evident that the fractionalized spin excitations 
are "hosted" in some form of solitons, this time, kink like 
ones. The topological aspect, especially the attachment 
of a high energy "string" to some naively fractionalized 
excitations, indicates a possible catastrophic consequence 
of having such an excitation. The issue of fractionaliza- 
tion therefore is intimately associated with the termina- 
tion of such a long range singularity, by dynamical means. 
This was also the idea behind the early work of Kivel- 
son, Sentha and Rokhsar, Rokhsar and Kivelson on the 
2D antiferromagnetic spin-1/2 systems in square lattices, 
and most recently Moessner and Sondhi on 2D frustrated 
antiferromagnetic spins in triangle lattices |lS| , pl] , ^7| ■ 

In the quantum dimer model in square lattices, again 
each dimer represents a singlet pair. To make contact 
with spin liquids of half spins, one also has to apply a 
hard core constraint so that only one link emanating from 
a site is occupied by a dimer. The Hamiltonian is 

H = -J(| =>< HI + h.c.) + V(\\\ >< HI + | =><= |) 

(3) 

where the first term characterizes the energy gain from 
the resonance between two parallel configurations of 
dimers (horizontal and vertical pairs) . The last term sta- 
bilizes each pair of the parallel dimers when V is negative 
but distablizes such a parallel configuration when V is a 
positive. The competition between them results in two 
distinct phases. 

When J«F, the respulsion between parallel dimers 
is dominating and a VBC phase which contains no par- 
allel dimers emerges as a result of the repulsion. When 
J> F, a column phase where dimers sit on the top of 
each other prevails. In between, it is believed that at 
J = V, an equal amplitude resonating valence bond liq- 
uid exists as a point where, crudely speaking, a column 
configuration and a VBC are degenerate. Each phase has 
fourfold ground state degeneracy as one can easily visu- 
alize and both are fully gapped and spin incompressible. 

In both phases, the creation of a monomer, or a spin- 
1/2 excitation provoks a disturbance of the crystal long 
range order, even with four degeneracy. The disturbance 
which is caused by a frustration in any dimension higher 
than ID, is long range by nature. It again can be rep- 
resented by a string which is emitted from a monomer 
or half spin in the bulk and ends only at the bound- 
ary or another half spin. We are coping precisely with 
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the same situation as in ID spin 1 chain. We therefore 
intend to conclude that spin-1/2 excitations, sometimes 
coded as spinons in the literature are not generic objects 
in this model and only exists at one point J = V when 
the ground state is an equal amplitude superposition of 
a number of configurations. Otherwise, only spin one 
excitations similar to gapped spin wave excitations can 
be observed though at each lattice site the elementary 
spin is one half. This surprising result that the excita- 
tions carry spins doubling the spins of underlying parti- 
cles, causes major frustrations over the early believe that 
spinons with half spins are elementary excitations in the 
system. 

However, as recently pointed out, this frustration can 
be partially removed when the 2D antiferromagnetic sys- 
tems are geometrically frustrated, say due to geometries 
such as triangle lattices. In that case, the authors of ]27[] 
were able to demonstrate that in a fraction of the param- 
eter space the long range order found in either VBC or 
a Column phase is absent in ground states. This relieves 
the frustration of creating a monomer in the bulk of the 
system and makes spinons possible. We will not pursue 
further here but refer to the original work. 

Existence of spinions in other 2D antiferromagnetic 
systems was also suggested in a few later works p2| , ^3[ . 
In p2| | , the authors showed the free spinon excitations in 
a crossed-chain model of a 2D spin liquid; in pgfl , the au- 
thors investigated the stability of a fractionalized phase 
in a kagome lattice and arrived at conclusions consistent 
with those in p7j|. 

The other powerful way to describe the fractionaliza- 
tion in magnetic systems is to introduce gauge fields to 
characterize various induced interactions. And this is 
also the language we are going to employ in the discus- 
sion of quasiparticles in this paper. We will give a peda- 
gogical summary of the development and refer to original 
works when the details are concerned. The power of the 
field theory in the study of quantum magnetism perhaps 
was first fully appreciated in a work on one dimensional 
spin chain [ B4J , In this pioneering work, a topological 
term was derived to distinguish the spin chain of inte- 
ger and half integer spins, especially the excitation spec- 
tra. The Haldane gap observed in ID integer spin chains 
but not in half integer spin chain is a manifest trophy 
of the field theory approach. However, a massive em- 
ployment of the gauge field based approach in the study 
of many-spin systems perhaps was triggered by a series 
of papers of Anderson et al. |]Jl(J on the mechanism 
of high temperature superconductivities, see for example 
- 39 1 . It was envisioned that an antiferromagnetic 



insulator under doping first becomes a resonating valence 
bond state(RVB), a collection of spin singlet pairs, which 
eventually is responsible for the superconductivity; fur- 
ther more, spinons likely form a fermi surface even in 
the charge insulating limit. The RVB theory indicates 
neutral spin-1/2 fermionic spinons and so the existence 



of these spinon excitations which is a key feature in the 
now widely accepted mechanism becomes the most cru- 
cial issue §,|0|@. 

However, an analysis of the spin disordered states us- 
ing the field theoretical approach reveals a rather frus- 
trating aspect of the problem. That is, when the long 
range Neel order disappears as the superconductivity 
emerges, the spin fluctuations effectively lead to U(l) 
gauge fields. And spin-1/2 but bosonic spinons always 
carry unit charges with respect to the induced U (1) fields, 
reflecting the Berry's phase a spinon develops under a 
spin rotation. So they interact with each other via U(l) 
gauge fields. The key feature of the induced C/(l) field is 
that in 2D even at zero temperature, they always mediate 
a confining force between two U(l) charges. Namely, the 
energy of two spinons separated at distance L is propor- 
tional to L and this catastrophe manifests itself by forbid- 
ding any excitation carrying a U(l) charge. Within the 
frame work of that approach, neither neutral fermionic 
spinons nor bosonic spin 1/2 spinons can exist. This 
is clearly reflected in the paper by Wiegmann jlSj ]. An 
interesting phenomcnological proposal to avoid this dif- 
ficulty is to introduce a Hopf term to have deconfining 
gauge fields and statistical transmutations at the same 
time j40|. The other interesting but more recent pro- 
posal was made by Senthil and Fisher; for discussions 
on this Z2 gauge field based approach, we refer to their 
original paper p6[. 
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Fig. 2 Following |20|(see also |62j) in disordered limits, a) 
Fourfold degenerate VBC states for "spin" S = 1 in a 
2D square lattice; b) Twofold degenerate VBC states for 
S = 2; c) Nondegenerate VBC states for S = 4; d) The 
creation of spin-1/2 excitations in a 2D VBC state. 

The Hopf-term-based proposal, promising and intrigu- 
ing at the first sight was proved to be inappropriate in a 
antiferromagnetic system where the spatial parity is un- 
broken |62l|63| . Attempts to derive the Hopf term micro- 
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scopically in Heisenburg antiferromagnetic systems |6^] 
indicate that instead there are other local Berry's phase 
terms which define the nature of disordered states. These 
Berry's phase terms result in destructive interferences be- 
tween space-time monopole events at different points at 
square lattices. For odd integer spins all monopole events 
with unit charges are suppressed, and only monopoles 
with double charges are allowed; for half integer spins 
both unit charge and double charge monopole events are 
forbidden and monopoles are always quadrupled. The 
ground state in this disordered limit has fourfold degen- 
eracy for half integer spins and twofold degeneracy for 
odd integer spins [|62| . This fascinating result was later 
substantialized and refined in [^o| . The authors of (2(]] in- 
troduced an SU {N) representation and explored system- 
atically the Berry's phase effect in the context of an ef- 
fective compact U(l) field theory. It was further pointed 
out that those degenerate states are of spin Peierls type 
breaking the lattice symmetry and correspond to the 2d 
valence bond crystal states closely connected with those 
early works of AKLT. They also emphasised that due to 
the confining property of the U(l) gauge fields certain 
topological orders are absent and the bosonic spinon ex- 
citations are confined. In another interesting work, these 
authors examined a frustrated Kagome lattice and sug- 
gested that the spin one half spinons are deconfined due 
to a charge 2e Higgs field disregarding the underlying 
spins p5[ . 

Though a Hopf term induced statistical transformation 
is believed not to take place in usual antiferromagnetic 
systems, because of the absence of time reversal symme- 
try breaking and parity breaking, theorists do agree and 
believe it should happen in systems such as chiral spin liq- 
uids [^2j. Then one can expect that the spin 1/2 spinons 
can exist in the excitation spectrum. As we should see 
below, a 2D p-wave superconductor is precisely such an 
example. The existence of a topological term in p-wave 
superconducting states was previously demonstrated in 
. Most attentions in those early works have been paid 
to many exotic properties of topological excitations in 
symmetry broken states. However, the issues of spin or- 
dering and spin-phase separation haven't been explored; 
in particular, spin fluctuations' effects on quasiparticles 
remain to be understood, especially under the influence 
of the Hopf term. This will be the focus of this article. 

In this paper, we will demonstrate the existence of 
spin rotation symmetry restored p-wave superconducting 
states and explore the properties of quasiparticles and 
collective excitations in these states. In section II, we 
scrutinize the spin ordering in the context of spin-charge 
separation, and prove the coexistence of a long range 
phase order and a short range spin order in a 2D spin 
triplet p-wave superconducting state. This spin rotation 
invariant state is characterized by a finite range spin cor- 
relation and hc/4e vortices as the elementary topological 
excitations. 



In section III, we characterize the topological orders 
in rotation invariant states. We show in a quantum 
spin disordered case the zero energy topological charges 
are conserved due to a quantum protectorate. In sec- 
tion IV, we examine the spin fluctuations' influences 
on excitations, particularly, quasiparticles and quan- 
tum collective spin excitations in 2D spin disordered p- 
wave superconductors. The elementary quasiparticles are 
Bogolubov-De Gennes (BdeG) quasiparticles hosted in 
gapless skyrmions and obey bosonic statistics in quan- 
tum disordered cases. The spin collective excitations 
are shown to be fractionalized ones carrying a half spin 
and obeying fermionic statistics, by contrast to the spin 
wave excitations in spin ordered p-wave superconduct- 
ing states (SOpSSa). In thermally disordered cases, BdeG 
quasi-particles form bound states. Spin triplet super- 
conducting states are believed to exist in 3 He, many 
heavy-fermion superconductors and most recently layer 
perovskite SY2-RUO4 crystals We believe the in- 

teresting properties studied in this work can be observed 
in certain limits. 

We should remark that the issue of fractionalization 
in spin triplet superconducting liquids was recently ad- 
dressed by Demler et al. J4l)] . There are two important 
differences between the current work and that of [ Ell . 
First, we are going to cope with an order parameter liv- 
ing in space [S 2 x S 1 ]/Z2 while Demler et al. dealt with 
[S 1 x S 1 ]/Z2- Though from the point of view of induced 
Z2 gauge fields which are the focuses of that work, this 
distinction is not important, it is however vital when spin 
textures and the Hopf terms are concerned. Second, the 
important aspect discussed in this paper, the Bogolubov 
quasi-particles hosted in skyrmions, depends crucially on 
the topology and is unique in the model studied here. 
The properties list in Tab.l. are novel ones of the spin 
disordered superconducting liquids. On the other hand, 
our theory is valid when the Z2 gauge fields are frozen 
and 7r-disclinations are gapped. For a detailed discussion 
about Z2 gauge fields and their effects on fractionaliza- 
tion, we should refer to the work by Demler et al. . 

Fractionalization, though has been discovered in very 
diversified strongly correlated condensed matter systems 
which bear no similarities at first sight, appears to share 
one remarkable common feature. It is the topological ex- 
citations interacting with electrons one way or the other 
which make exotic quasiparticles or collective excitations 
possible. This also lies in the heart of earlier examples 
discovered in field theoretical models and mathematical 
physics PJ^,p|,^2t . Indeed, more sophisticated considera- 
tions have led us to believe that quantum number of ex- 
otic excitations seen in quantum Hall effects or appearing 
in theories on spin liquids should depends on the ground 
state degeneracy of the liquids in certain multiply con- 
nected manifolds. The connection between the ground 
state degeneracy and the fractional charges in quantum 
Hall liquids was raised in p4 22 1; the topological orders 
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in spin liquids were emphasised in a few occasions, see 
for instance p3| , p8|] . As we will see, this notion plays a 
paramount role in our case as well. 



II. SPIN-PHASE SEPARATION AND SPIN 
CORRELATIONS 



which satisfies the usual commutation relationship 
[S Q ,S^] = ihe a p 7 S-y. Spin S states are expressed as de- 
gree 2s polynomials of u and v; u 2S ~ m v m , m — 0, ...2S. 
The scalar productor is defined as J f*(u,v)g{u,v)d£t. 

As an exercise, one calculates the kinetic energy 
E(fl, X ) = (87T)- 1 J dtt u*(ft)^-u(fl) of the spin one 
half state u(fl) as 



A. Spin-phase separation 

The term "spin-phase" separation is used in p-wave 
superconducting liquids to represent a separation of the 
phase dynamics and spin dynamics. So the low lying 
goldstonc modes can be well classified as phasons and 
spin waves. However, this does not imply, at least not 
directly, the coexistence of short range spin correlation 
and a long range order in phases. In contrast, in the mean 
field approximation of BCS type, the "spin-phase" sep- 
aration defined above always indicates a spin long range 
order if a phase rigidity is present. So the fact that we 
can have spin disordered superconducting liquids in 2D 
should not be taken for grant. And indeed, in any three 
dimension sample both the spin order-disorder transition 
and superconductor-metal transition take place at tem- 
peratures within the Ginsburg region and practically si- 
multaneously. 

To get oriented with these statements, let us first take 
a phcnomcnonlogical approach. We are going to consider 
the cooper pairs condensed in spin triplet states, or spin 
one states. The spin component of the pair wavefunction 
can be written in Wely-Majorana representations as 



«(ni)«(n 2 ). 



(4) 



Here u(fl) and v(fl) are expressed in terms of a vector 
(it, v) on a unit sphere, 



E(tl,x) = — |V X 
2 m 

A = cos0V^,F M!/ 



A| 2 
1 



2m 1 



ivnj- 

dflg $ri„ 



, $7 . """ft t/J " 7 

2 a Pi a q x q Xv 



(7) 



where F^ = e^d^A^. 

Therefore, the charge degree of freedom \ is coupled 
to a topological gauge field F M1/ . It is this coupling that 
in general fails to yield a spin-charge separation in anti- 
ferromagnetic spin systems. Eq.[?] also indicates that the 
current carried by a state u is 



j(n) = i(vx-cos0^) 

m I 



(8) 



with the last term from a spin rotation. In terms of the 
current, a spin-phase coupling manifests in an extra con- 
tribution due to the spatially varying n. 

The supercurrent J s (J7i, r2 2 ) carried by a paired state 
can be obtained by calculating 

j,(n 1 ,n a ) = 6 / ^u*(ni)»*(n 2 )^tt(niMna). (9) 

J 47r Zm 

However, a direct calculation shows that J = if fix = 
0,2 as in a quantum spin nematic state. Otherwise, 



j.(ni,n 2 ) 

with 



1 



1 



Vx + — sin 6[66S/<f> ~ 5<j)V9] (10) 



u(fl) = exp(zx)[cos - exp(-z-)u + sin - exp(i-)v ] 
v(fl) = exp(-ix)[-sin^exp(-i^)u + cos ^ exp(i^)w ] 



(5) 
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2 

h+&2 



2 

h + 4>2 



(11) 



where fi = (sin cos 0, sin sin 0, cos 0). And the unit 
vector (u,v) can be written as u = cos#o/2exp(i0o/2), 
v = cos#o/2exp(— icj)o/2). The wavefunctions in Eq.4 
form an overcomplete set fro spin-one states. The quan- 
tum spin nematic state which interests us below corre- 
sponds to a state S~2i = fi 2 . But to see how phase-spin 
separation takes place we keep two unit vectors indepen- 
dent at the moment. 

In this representation, spin operator is defined as 



S+ 

s z = 



■■ u—, = v — , 
ov ou 

1 d 8 
~[u- v—\ 

2 ou ov 



when 56, 6(f) < 1. Or 

J s (n+ ^n- ~) = -^-eap^naXpV^. 



(12) 



(6) 



From here one concludes that there is a coupling be- 
tween the spin and charge sector. However, the current 
described in Eq.[f| is much smaller than that in Eq.|| as 
1 approaches zero. As a result, the effective gauge field 
in this case is 

F^(n,l)^- eQ , 7 [n Q ^^ + — 1^]. (13) 

It is negligible compared with the usual Pontryagin field 
in Eq.^ because of its dependence on 1(<C 1). 



G 



Indeed, a direct calculation of the kinetic energy of the 
cooper pairs condensed at a u(ft + l/2)v(£l — 1/2) state 
shows 

E = ^|Vtt| 2 + ^|V X | 2 + • A*(n, 1) (14) 

2m 2m m 

and A p is the vector potential of F p which vanishes as 1 
approaches zero. p s is the superfluid density. In the low- 
est order gradient expansion, we conclude a spin-phase 
separation as far as the linear dynamics is concerned. 

Having said this, we should also emphasis that the 
pairs which can be considered as spin one particles carry 
no U(l) charges denned with respect to the Pontrya- 
gin field F^,,, because of the spin-phase separation. As 
we will see later, on the other hand, the Bogolubov- De 
Gennes quasiparticles do carry unit U(l) charges so that 
spin-charge separation does not take place. 

Eq.|l4] further indicates that in the mean field theory, 
the phase rigidity always induces a spin rigidity and the 
superconductor- metal transition is always accompanied 
by an establishment of a magnetic order in n. A more 
sophisticated approach should take into account the fluc- 
tuations in both spin and phase sectors. At a dimension 
higher than two, one can easily confirm that 

1 ■ (15) 



pom 3 

where T c sm is the superconductor-metal phase transition 
temperature, T* is the induced spin order-disorder tran- 
sition temperature and T c is the mean field critical tem- 
perature, po is the total density and £(0) is the zero 
temperature coherence length. Eq.[l5] shows the sub- 
tlety of the spin-phase separation and consequently the 
coexistence of short range spin correlation with a long 
range phase order in a 3D system. The righthand side of 
the equation which virtually characterizes the size of the 
Ginsburge region is always much less than one. There- 
fore, in most of spin triplet superconducting states dis- 
covered so far, the mean field picture is more or less cor- 
rect and from a practical point of view, the two tran- 
sitions do happen at the same temperatures in a bulk 
system. 

However, in general, the spin disorder and the phase 
disorder are driven by very different mechanisms and 
the two phase transitions belong to different universality 
classes. We should show that Eq.^H] is false when a 2D 
p-wave superconducting liquid is concerned and a short 
range spin correlation does coexist with the superconduc- 
tivity. For that purpose, we now turn to a microscopic 
derivation to include both the zero point motions of n, x 
and spectral flows in fluctuations. 



B. The effective action 

We will be specifically interested in a p-wave supercon- 
ductor with an order parameter 



A(k)=A (T)(o5cr.n)(A x + ilfc tf )eK P (*x), (16) 

where I = ±1, and A is defined in the space of spin | ,| 
^2| , ^3| . n is a unit vector living in a two-sphere, S 2 . 
This superconducting state breaks a U(l) symmetry of 
the phase of the order parameter, a U(l) symmetry of 
the rotation along z axis and a S 2 symmetry of the ro- 
tation of spin quantization axises. In addition, it breaks 
the pairty invariance connected with I = ±1 and a time 
reversal symmetry. 

One can always undo a rotation along the z axis by 
a £7(1) gauge transformation and restoring this rotation 
symmetry also means an absence of superconductivities. 
We will not be interested in this symmetry restoring but 
focus on the rotation symmetry breaking or restoring in 
the spin sector within the superconducting phase. For 
this reason, we reserve the "rotation symmetry" solely 
for the S^-symmetry in the spin sector. Consequently, 
the rotation invariant state suggested in the title refers 
to a state where the S 2 spin rotation symmetry is re- 
stored. 

The corresponding Hamiltonian for the superconduct- 
ing state in the Nambu space of '5 = (tp + ,iT2tp) can be 
written as, 



H = a 3 e+ °i{ d i>&}- 



(17) 



where A is defined as A = Ao exp(i(73x)(n • t) and 
e(k) = h 2 \i 2 /2m — €f- We use a as the Nambu space 
Pauli matrix and r as the spin space ones. We assume 
the spin-orbit interactions are weak and n is a unit vec- 
tor in a sphere S 2 . The internal space of the symme- 
try broken state is 72. = [5 1 x S 2 ]/Z2, similar to that 
of Bose-Einstein condensates of 23 N a studied re- 
cently. The order parameter observes a discrete symme- 
try: A(n, x) — * A(— n, x + 7r) and represents a quantum 
spin nematic p-wave superconducting state. 

To obtain an effective theory, we integrate over the 
fermionic degrees of freedom and make a gradient ex- 
pansion. At low temperatures, we obtain the following 
results 

C = 5 ab P so P s6 - T $ 2 + -^(E 2 + B 2 ) 
+«S^V a n a V i ,n /3 - T a "d t n a d t n p + — e^W 3 ^. (18) 



and <I> = IdtX + e<f> e 



In Eq|J, P s = iV X + eA e 
are the gauge invariant momentum and potential respec- 
tively. The superscript em is introduced to distinguish 
the usual electric magnetic vector potential A em from 
the topological field A defined in terms of n in Eq.0. 

And W 3 is the vector potential of F^ (defined in 
EqM). 
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(19) 



In deriving Eqs.|18|,[19|, we have introduced two vector po- 
tentials denned in terms of a U(l) rotation U c and a 
SU{2) rotation U s , 



A CM = iU c 1 u c 1 d ii uJJ c = cr 3 {Al m + i<9 M x), 



(20) 



(/j = 0,1,2, stands for coordinates in 1 + 2 dimension 



space.). u c — exp(— i(i3^&-), and u s — exp(iT3^p). 

The time dependent gauge rotations a(t),f3(t) are in- 
troduced to preserve the anti-periodic boundary condi- 
tions for Fermions under a nontrivial gauge transforma- 
tion p9j. And the two gauge rotations are defined as 



Pit)- 



U+U s = 1, U+d ■ tU s = T 3 , 
U+U c = 1, CZ+fTi exp(icr 3 x)£4 = cr 4 



(21) 



At last, the temperature dependences of all coefficients 
are given in Appendix B. Let us define S^ = 5 al3 Sab, 
and dimensionaless coefficients 



Po 5 



T) — VqTq, S a b = S a b, T a b = 5 aP VqT . 

2m 



(22) 



Vq is the averaged density of states at the fermi surface of 
the normal state, po is introduced as the density; m and 
€f are the mass and the fermi energy respectively. Then 
7q, S a b, T and M are calculated as 



To 



1 



d 2 k 



2kT 



4, J ( 2 7r ) 2 - F2,k,, ' ,!l( ^ , • 
2m /" d 2 fc 

( 

d 2 /fc 



Sab = / 7^2V V6Fl(k), 

Po J (2tit 



voJ (2tt) ; 



^i(k); 



4tt £ 3 (k) 5^ dk t 



Fi(k). 



(23) 



all, the dynamics of spin n and phase x is completely 
separated at the low frequency limit (except the entan- 
glement due to the Z-i group in the functional inte- 
gral [ ^H|4^j52] , |53]] , which we will not discuss in this pa- 
per.). It reflects a spin -phase separation in the p- 
wave superconductors, as argued in the previous sec- 
tion on a phcnomenonlogical ground. Second, for an 
isotropic fermi surface which interests us in this arti- 
cle, Sab — S a bPs{T)/2m, and p s {T) is the temperature 
dependent superfluid density which vanishes at the crit- 
ical temperature T c . So the spin and phase dynamics 
are characterized by an 0(3) cr-model (NLctM) and an 
xy model respectively. At the mean field approximation, 
n = e 2 and x = const.. This corresponds to a conven- 
tional SOpSS. There are three Goldstone modes; two of 
them are spin waves = (1, ±i, 0) with a linear dispersion. 
In an isotropic case, S a b — 5 a bS; the spinwave velocity 

is v s (T) — vp\l S~/2nT. And the last mode is the usual 



plasma wave, with a dispersion u) = yj2ne 2 pok/m in 2D 
at T = 0. 

The topological term was previously derived in [Q|35| . 
There, J\f is effectively the number of skyrmions living in 
the (k x ,k y ) plane if we define a skyrmion configuration 
in the plane as M(k) = (vAkr, i>Ak y , e(fc)) . Namely, 

* = §s F ^E-\\<i)e a ^M. a d k Mpd ky M. 1 and N = 1 
for a single band case. This term determines the topo- 
logical orders in the fields F^ u and defines the structure 
of quasiparticles. We investigate the spin ordering, quasi- 
particles and topological excitations based on Eq.fl8[ 

The topological term arises naturally in a system where 
the parity and time reversal symmetry are broken. It is 
therefore not suprising similar terms can be present in 
the effective theories of other unconventional supercon- 
ductors where these symmetries are broken. Some im- 
plications of this phenomenon on (spin) Hall effects were 
indeed explored recently in a series of papers |36] ^9) . We 
refer to these original papers for a detailed discussion on 
transports. 



E{k) — y/e 2 (k) + kZAy/kp, and A is the temperature 
dependent gap. The director M is defined as M(k) — 
(va^x, VA^y, e(fc)) ; va = Ao/kp- And F\p are calcu- 
lated in appendix B. 

Following Eq.^3], TqVq 1 is a unity at T = and varies 
smoothly as a function of the temperature |48| , Vq is the 
averaged density of states at the fermi surface. S a b on the 
other hand approaches zero at the mean field critical tem- 
perature T c . When the director n is a planar vector con- 
fined in the = tt/2 equator, i.e. n = (cos </>, sin </>, 0) , the 
spatial gradient terms coincide with previous results for 
the A phase of 3 He [^0|,^l). The action is valid when the 
frequency and the wave vector are smaller than A (T) 
and £n_l( T ) = A (T)/v F respectively. 

Eq.nS suggests a few important properties of the quan- 
tum spin nematic p-wave superconductors. First of 



III. SPIN DISORDERED SUPERCONDUCTING 
LIQUIDS 

Because of an extra branch of Goldstone modes in the 
spin sector, the spin order is more fragile than the phase 
order in the problem when fluctuations are taken into 
account. This provides a possibility of having spin dis- 
ordered p-wave superconducting states (SDpSSs) where 
the S 2 symmetry is restored and only the U(l) symmetry 
is broken. Such a state which is spin rotation invariant 
in nature differs from the conventional SOpSSs where the 
S 2 symmetry is broken and there is a long range order 
in n. 

The finite temperature phase diagrams of the 0(3) 
Nonlinear a Model (NLcrM) were previously analyzed 
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in great details |3(J. In the current situation, just as 
the superfluid velocity p s (T), all coefficients in the ac- 
tion, <5>"jf , T a b, S a b, T depend on temperatures because of 
quasiparticle excitations. Taking this into account, we 
arrive at the following results in 2D. 

QSDpSS When A 3> £_f, the spin order is spoiled 
by quantum fluctuations and there is a short range spin 
correlation even at zero temperature with a finite corre- 
lation length. 1 This state will be named as a quantum 
SDpSS, or QSDpSS. 

TSDpSS When Ao -C ep, the spin order is established 
at zero temperature and the correlation length 



A S (T) 



,A S = Texpi 



27r[p s (T)/2m-r] 
T 



(24) 



is finite only at finite temperatures (in a saddle point ap- 
proximation) . Here T ~ Ao (T) . This state is a thermally 
SDpSS, or TSDpSS. Discussions about quantum orders, 
zero energy skyrmions and BdeG quasi-particles will be 
carried out paralelly for both cases. 



A. Instantons and their suppressors 

The nature of spin fluctuations in a QSDpSS can be 
explored in a spinor representation of Eq.[l8| By intro- 
ducing r] + T7] = n, 77 — (771, t/2) t and ?/ + r; = 1, we obtain 
for 77 the following Lagrangian in SDpSSs, 



A )n\ 2 I As ^ T K + n I — f^ a A F i 

(25) 



And 77 is a bosonic field carrying a unit charge with 
respect to A fields and spin 1/2. In Eq.||[ 2/ 2 = 

2mAo / V ST pq ; we have introduced the following rescal- 
ing: t —> it£o/v s , r — > r£o and Eq.^5| is written in a 
(2 + 1) Euclidean space. 

When the spins are disordered, the gapped spin fluctu- 
ations induce an effective Maxwell term. Upon integrat- 
ing over spin wave excitations, based on a general consid- 
eration of the gauge invariance and the parity symmetry 
breaking, we conclude the NLctM should be reduced to 



1 

~9 



4tt 



W J F 



(26) 



<?(A S ) is a function of A s the spin gap. A direct dia- 
grammatic calculation yields the following estimate for 
9 



Ap 
8ttA, 



(27) 



In (2 + 1) space x = (r, r), it is convenient to introduce 
a field, H, ( = ^'""T^. H T = F xy represents C/(l) mag- 
netic fields along z direction, H x = F yT and H y = F T:C 
are the y, x-components of the electric fields. H is the 
solution of the Gauss equation 



d v H v (x) = 2J Q m S(x 



(28) 



l)d 



in the presence of space-time monopoles {x m } in (2- 
Euclidean space. 

Following Eqs.^6], |28| one can obtain an effective ac- 
tion in terms of monopoles' coordinates. However, the 
Hopf term leads to destructive interferences between ro- 
tated skyrmions terminated at space-time monopoles. 
All monopole events are therefore completely suppressed. 
This is reflected in the action for monopolc-like instan- 
tons. But before presenting the result, let us characterize 
the topological structures of these instantons. 

Consider a rotating skyrmion terminated at the origin 



S space 



n(p, 4>;t) 



where 



0(p,r) 



sin a cos 0, sm t) sin <p, cos ( 

6-7(7-), 

larcos — = (t) , 



(29) 



Vp^T 



1 (t + L t )-j(L t ) = N2tt. 



(30) 



To obtain a desired boundary condition for 77-quanta and 
quasi-particles discussed in the next section, we also re- 
quire 7(t + L T ) — 7(t) = N2tt. The perimeter along r 
direction L T is finite for the discussion of a finite temper- 
ature case; L T = T _1 . 




x This is supported by a mean-filed calculation in a CP N 
representation. 
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Fig. 3 Skyrmions terminated at space-time monopoles at 
different points carry different phases 7b, which result in 
destructive interferences. 

The gauge fields in this "texture" are 

= j-cos9(p), A T = i cos6»<9 t 7(t); 



H r = 0(r) 



2Qi 



(p 2 



r 2\2- 



H p = 6(r) 



(/9 2 +T 2)2 

2Q P t 



S(t) 



Qp 2 -r 2 



' (/9 2 +r 2)2 

One can easily confirm that 

V H = Q2ir5{r)6(r), 



2pp 2 



(31) 



(32) 



representing a monopole of charge Q. So a space-time 
monopole event terminates a skyrmion. As p, r approach 
infinity, H(p, t) becomes vanishly small. The action of 
this Euclidean space monopole charge is finite(a ~ 1), 



Sn 



j- J drd 2 rU 2 (p,r) 



aQ 2 _ 



(33) 



(a is a constant of unity.) in addition, it has a a phase 
factor following the Hopf term, which characterizes a ro- 
tation of the skyrmion during its duration 



7b(t ) 



QAf 



i{ l t) - 7(t"o)]- 



(34) 



Note that the term inside the bracket is not subject to 
the periodic condition imposed for the values of 7 in 
Eq.[30|. 7s obviously depends on the temporal coordi- 
nate at which the skyrmion is terminated, introduced 
explicitly as tq in Eq.|34|. 

Taking into account this feature, the action of space- 
time monopoles {Q a , N a } centered at r Q , r a can be writ- 
ten as in Eq.Cl. The most important feature thus is the 
space-time monopoles positioned at different points carry 
different phases 7s(Q q ,t q ,). It is interesting that for 
Heisenburg antiferromagnetic spins at square lattices, the 
space-time monopoles carry spatially dependent Berry's 
phases which lead to spin -Peierls ground states |]2| . 

A simple calculation in a saddle point approximation 
shows the energy density 



E 
1? 



-<5(A0^ex pi 
So 



167T<7 



cos(x s ). 



(35) 



Eq.^35| indicates that only when Af = 0, monopole events 
are liberated and the ground state is \s — n2n, n = 
0, 1, 2.... However, for Af ^ the ground state is of arbi- 
trary Xs with infinite-fold degeneracy. The suppression 



of space-time monopoles is a unique feature of QSDpSS 
which we are going to explore in some details. 

To appreciate the robustness of Eq]3|, we offer an al- 
ternative view from the perspective of the spectral flow. 
The absence of space-time monopoles in QSDpSSs is a 
direct consequency of a fermionic zero mode. Support 
for this conclusion comes from the following argument 
similar to that of the Witten's global anomaly |l3| . The 
Hopf term in the action, 

Af r , . 1 



hopf 



4tt 2 



(A M + -d t p6„ t0 )d u A x , (36) 



suggests that the determinant of the fermionic operator 

C = d T -T-W -H(Wi), (37) 

transforms nontrivially under gauge transformations. Let 
us classify a gauge transformation as 



W,(r,r)^W (1 (r,r) + (/„+9,[/ n , 



(38) 



where U n belongs to a nontrivial element of ^(S 2 ) group, 
or a Hopf texture. This has to be supplied with 



(3(L T ) - P(t + L T ) = n2ir. 



(39) 



Under the transformation defined in Eq.[38[ the determi- 
nant acquires a negative factor due to the Hopf term 

(-l) Mn det£. (40) 



detC 



To see how this can happen, we analyze the spectral 
flow during different stages of gauge transformations. We 
introduce an additional parameter a to monitor the spec- 
trum of the Lagrangian operator such that 



W M (r,- 



0, a 
W»(r,r), a 



(41) 





(a) (b) (c) 



(d) 



(e) 
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Fig. 4 The different stages of gauge transformations in 
Eq.pl]. a) a = 0, b)| > a > 0, c) a — |, d)l > a > |, 
and e) a = 1. A skyrmion (indicated by the tubes) is 
terminated at a monopolc-antimonopole pair (indicated 
as filled circles). In f), a level crossing in the quasiparti- 
cle spectrum takes place at a = 1/2. 

At a — 0,1, the transformation is a pure gauge and 
the spectrum is given as 



(42) 



oj n L T = 7r(n + l/2). Recall E(k) = ^ A 2 + e(k) 2 and the 
spectrum of the Lagrangian operator is fully gapped. 

In between, we have situations shown in Fig. 3. We 
start with a skyrmion which is emitted from a monopolc 
and terminated at an anti-monopole. Slowly we send the 
two monopoles to r = ±L T /2. At a = 1/2, we will have 
two monopoles inserted into the boundaries of 2 + 1 D. 
The gauge field W at a = 1/2 precisely corresponds to 
that of two space-time monopoles. As two monopoles 
approach each other, the skyrmions are rotated by 360° 
degrees. 

To obtain a minus sign as indicated by the Hopf term, 
the spectrum has to cross each other once when the pro- 
cedure specified in Eq|3^ is carried out. Since the gauge 
fields are symmetric with respect to a = 1/2, the spec- 
trum is symmetric as well. The crossing has to take place 
at a = 1/2 to give us a designated —1 contribution as 
suggested by the Hopf term. 

This is the point where the gauge fields are most sin- 
gular, defined by the monopole fields and the spectrum 
of the Lagrangian operator can be closed. The probabil- 
ity to have space-time monopoles is proportional to the 
ratio between the fermionic determinant in the presence 
of the monopole detCi and that in a trivial background 
detCa, 



detCi 
detCo 



IL 



£ 2 (k))' 



(43) 



Since the eigen value of the operator Co is fully gapped, 
the zero energy mode flo = in the presence of a space- 
time monopole completely shuts down the monopole pro- 
liferation. This completes our argument. The argument 
only depends on the general principle of gauge invariance 
and not on the details of the gauge transformation. 

It is important to mention that the action Eq.|2^ is valid 
when the time scale in the problem is shorter than A" 1 . 
The duration of the monopole should be understood as 
A s . The action for an individual space-time monopole 
in Eqs.[53],[54] is well-defined only when the perimeter L T 
is longer than Aj 1 ; otherwise, the thermal nucleation 
dominates. The crossover happens when 



T ~ T* 



1 



A S (T)' 



(44) 



It is therefore clear that in the SDpSS driven by quan- 
tum fluctuations, QSDpSS, the suppression of monopoles 
is extended into a finite temperature. On the other hand, 
for the SDpSS driven by thermal fluctuations, TSDpSS, 
T is always higher than T*. The Skyrmions are ther- 
mally nucleated; the topological phases are irrelevant for 
the discussions and the suppression is absent. The gauge 
fields and excitations in QSDpSSs are qualitatively dif- 
ferent from the those in TSDpSSs. 



B. Zero energy skyrmions and quantum orders 

To understand the issue of the conservation of 
skyrmion charges in SDpSSs, we first exploit at least 
three important intraconnected consequencies of the zero 
modes. These three important features are vital for the 
discussions of quasiparticles in a SDpSS. First, the fluc- 
tuations of space-time monopole events per unit volumn 
defined as < Q 2 >= L~ 2 d 2 E(xs)/dx 2 s are 



<Q 2 >= W ^exp(--^- 



)• 



(45) 



in QSDpSSs. Eq.^ shows that at any finite J\f, all 
monopole events are suppressed due to destructive inter- 
ferences between monopoles with different rotation an- 
gles 7b. In TSDpSSs for the reasons mentioned at the 
end of the last subsection, the suppression is absent. 

Second, Eq.|45| indicates the conservation of the 
Skyrmion winding number at Af ^ in a QSDpSS, that 
is in the absence of a spin rigidity. To see this, one first 
notices that the equation of motion for is 



9Hfc P* a i a 



(46) 



Eq.[46| demonstrates that the topological charge is con- 
served if singular space-time events are not allowed so 
that the product n ■ dn in the right hand side of the equa- 
tion vanishes identically. And the topological charges 
fluctuate only when singular space-time events are per- 
mitted in the ground state. 

This can be more explicit in 2D. We define 



CV/(i)({n(r)}) = -L J dxdyU Z7 



(47) 



as the total number of Skyrmions living on the 2D sheet. 
As shown before, a space-time monopole essentially con- 
nects a trivial vacuum to a Skyrmion configuration. In- 
deed, following Eqs.|^,||, 



dC uw (T) 
dr 



Q m 5{r - r m ), 



(48) 
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where the surface contribution has been neglected since 
we are interested in the leading contribution to Cum 
per unit square. Following Eq.^8|, each monopole event 
causes a change in the topological charge Cum by one 
unit. The skyrmion charge is conserved in the absence of 
the space-time monopoles in (2 + 1)D. 

Combining the results in Eqs.[l5[ we conclude that 
the skyrmion charge is indeed conserved in QSDpSSs. A 
skyrmion is a well-defined topological configuration even 
in the absence of spin stiffness, a remarkable result which 
only holds at M ^ 0. Eq.^3] indicates that the zero mode 
in instantons serves as a quantum protectorate to pre- 
serve the topological charge conservation. The energy 
of skyrmions originates from their interactions with spin 
fluctuations. For a size A skyrmion, following Eq.^6| we 
arrive at an estimate 



E sk . = ^Ao(T) 



(49) 



which vanishes as A goes to infinity. 

Thirdly, this interesting destructive interference be- 
tween different monopole events leads to very distinct 
behaviors of fields F M1/ . Formally, one can introduce an 
order parameter of the Wilson loop type to characterize 
the induced U (1) fields. The Wilson-loop integral defined 
as 



W u{1) = (Vexp(il A-dr)) 



(50) 



has different asymtotical behaviors in the large loop limit 
in the presence or absence of the topological order in 
Cum- When the topological charge Cum is conserved, 
the monopole-like instantons which connect topological 
different configurations are forbidden; the exponent in 
the Wilson loop integral is a linear function of the perime- 
ter L of the loop, i.e. 



W u{1) = exp(-iCi) 



(51) 



and the gauge fields are deconfining. Eqs.^5 48 shows 
that in QSDpSSs, the gauge field is deconfining. 

In TSDpSSs, on the other hand, Cum is not conserved; 
the instantons are allowed and one can confirm the ex- 
ponent of the Wilson loop integral is a linear function of 
the area S enclosed by the loop, i.e. 



Wum = exp(-SC 2 ). 



(52) 



and the gauge fields are confining. For most of layered 
p-wave superconductors, given the gap energy at zero 
temperatures Ao ~ IK and the Fermi energy of order 
leV, we conclude the superconductors are spin ordered 
at the zero temperature but could be in TSDpSSs (in the 
absence of spin-orbit coupling) as shown in Eq.^. TS- 
DpSSs can be studied in many layered structure super- 
conductors. We will be concerned with both situations, 
the TSDpSSs and QSDpSSs in the following discussion. 



IV. QUASI-PARTICLES IN A SDPSS 

A. U(l) charges of Quasiparticles 

We now employ the generalized Bogolubov-De Gennes 
equation to study the properties of quasiparticles in a 
SDpSS. In the presence of a topological configuration of 
n(r), it is convenient to introduce a gauge transformation 

- exp(^a 3 ) exp(z^T 3 )[/ s (n)(7 c (x)* (53) 

and work in a rotated representation; then one obtains a 
new Hamiltonian 



H = cr 3 e(iV - A c - A s ) + 

VA {(?iT3,idi - Ad - A si }- 
i=l,2 



(54) 



Here ua — Ao/kp. The phase factors of a, f3 are to 
ensure antiperiodic temporal boundary conditions for 
quasi-particles, \P(t) = — ^(r + L T ). L T is the perime- 
ter along the direction of r; a(r) — a(r + L T ) = n2n, 
/3(r) - (3(t + L T ) = n2n. 
Explicitly, we can have 



Us 



cos 2 exp ( — i- 
sin | exp (ij 



■exp (if) 



(55) 



A direct calculation shows that 



w^a = iu+(e,^)d fl u s (e, ( t>) 

W 3 = I C os^ + ia M a. 



(56) 



It is clear that the W 3 field represents the vector poten- 
tial of F^, the Pontryagin field introduced before, 



W 3 = A, + ^ T /3(r)^, Q . 



(57) 



At last, the corresponding Lagrangian density is 

C B deG = * + (9r - u 3 A^ - T ■ A s0 - H( A™ , A sp )) * . (58) 

Therefore, the interaction between a quasiparicle and 
spin fluctuations can be best characterized in terms of 
gauge fields introduced above. Following Eq.[54|, besides 
a usual electric magnetic charge defined with respect to 
A e .m. fields, a BdeG quasiparticle also carries a unit 
£/(l)-charge with respect to A M fields and is minimally 
coupled with F M „. So both phase fluctuations and spin 
fluctuations affect the quasiparticles. 
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However, only the gradient of the phase, or superfluid 
velocity v s enters Eq.[54|. Though phase fluctuations are 
enormous in 2D, the fluctuations in v s 



(63) 



< v s v s >= V, 



2 T PO 

e F p s {T) ' 



(59) 



are small, except within the narrow Ginsburge region of 
the superconductor-metal transitions. We therefore con- 
clude the BCS spectrum is not affected by phase fluctu- 
ations at temperatures far below T c . 

The spin fluctuations in SDpSSs on the other hand de- 
fine all important properties of quasiparticles. A spin-^ 
excitation interacts with the spin fluctuations via gauge 
fields 



H i = \ J drJ ^ A M 

J M = eS(r - r a )^, J = eS(r - r a ). 

The energy of a spin-i excitation therefore is 



(60) 



E=~hi(Texp(-^J T ff/(f)rff)) (61) 
In terms of the induced connection fields, 
£ = --^-ln(exp(- ~ f rfx f dx 2 F 02 )) {FV} , (62) 



and express the topological term in a form of minimal 
coupling j^A^. By minimizing the action of L + CsdeG 
with respect to Ao, A em and taking J\f = 1 at a low 
temperature limit, we do obtain a saddle point equation 
4-7T < ^ + t 3 ^ >= e z V x A, < * + cr 3 5' >= 0. This shows 
that a skyrmion configuration carries a half spin but no 
charge. In other words, a spin i but chargeless BdeG 
quasiparticle is hosted by, or confined with a skyrmion, 
with the confinement mediated by the spin fluctuations. 
As a result, quasiparticles which are charge neutral with 
respect to U(l) fields can be created. It is the vanishing 
energy of skyrmions in spin disordered superconducting 
liquids which makes these topological excitations rele- 
vant to the discussions of the low energy quasiparticles. 
In spin ordered liquids, the skyrmions' energy can be as 
high as the fermion energy and skyrmions are irrelevant 
for low energy quasiparticles. 

For the study of a BdeG particle confined within a 
skyrmion in QSDpSSs, consider a skyrmion with nw — 1- 
In polar coordinates (p, (f>) , the director has a spatial dis- 
tribution as 

n(p, (f) — (sin 9(p) sin </>, sin6*(p) cos </>, cos 9{p)) (64) 

and 9{p) is an arbitrary function of p, the asymptotics 
of which is 9{p = 0) = and 9(p -> oo) = tt. The 
corresponding A field can be chosen as 



and the average is taken over the partition function 
ZdF^}). Only the y-component of the electric field 
defined by F 2 contributes to the energy. Say it differ- 
ently, the energy of a BdeG particle is determined by the 
Wilson loop integral of A fields. The evaluation of Eq.|62] 
depends on the hidden topological long range order con- 
sidered in recent preprints |53| (see also p8| for a general 
description on quantum orders), and provides all sorts of 
information on excitations. We will use this quantity to 
classify distinct structures of quasiparticles in TSDpSSs 
and QSDpSSs. 



B. The statistics of quasiparticles 

In QSDpSSs, all monopole configurations are sup- 
pressed because of the destructive interference caused by 
Hopf terms and < Q 2 >= 0. The Wilson loop integral 
decays exponentially as a function of the perimeter of the 
loop. The interactions between BdeG quasiparticles me- 
diated by the topological fields are perturbative and 
the BdeG quasiparticle energy is finite. But most impor- 
tantly, in this case, the skyrmion winding number Cfj^) 
is a conserved quantity; skyrmions themselves carry J7(l) 
charges with respect to the fields of A M . This is indicated 
in Eq.^ if we introduce the skyrmion density-current 
density as 



A = 1 - cos 6(p) c 



2psin0 " " " 2p dp 

To facilitate the calculation of the fermionic spectrum 
hosted in a skyrmion, we further assume that 9{p) = 
7t6(p - pa) is a step function of p. The SU(2) field W Q 
thus takes a simple form; 



(65) 



W?=A < (i = l ) 2),Wg=0; 
= W?, = 0. 



(66) 



The gauge field A M induced in the spin rotation above 
is simply the connection field of the Berry's phase. That 
is, J c drA = Q.c/2 where fic is the solid angle spanned 
by a closed loop C swept by the director n on the unit 
sphere, when the particle moves around a skyrmion. In 
fact, Eq.|65| represents a monopole field with the dirac 
string pointing at — e z direction. 
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Fig. 5 The composite structure of quasiparticles in QS- 
DpSSs. A bare BdeG quasiparticle (upward lines) and a 
skyrmion (doward lines) see each other as a dirac string 
pointing at =F e z, following Eq.|65|. 

When a spin- 1/2 BdeG particle with spin up in the 
basis J7 s (n)(l, 0) T moves in a closed circle in a skyrmion 
defect, the acquired Berry's phase of 7r[l — cosO(p)] which 
approaches 2n at an infinity p. Similarly, a skyrmion 
moving around a spin-1/2 particle acquires a Berry's 
phase of ir cos 6(p). The total phase involved in mov- 
ing the composite around the other one is 2ir. Conse- 
quently, under the interchange of coordinates, the two- 
body wavefunction of composite quasiparticles acquires 
an additional it phase because of hosting skyrmions, 

*(ri,r 2 ) = tf(r 2 ,n). (67) 

Eq.^ also follows the linking number theorem for 
skyrmions jPjJ. These composite quasi-particles are 
therefore Bosons. We also observe at e = 0, BdeG quasi- 
particles are charge neutral with respect to an em field; 
they also carry zero U{1) charges so to minimize the in- 
teraction between composite excitations. 

The life time of the composited quasiparticles is deter- 
mined by the minima of the life time of skyrmions which 
is limited by the quantum nucleation rate, and the usual 
inelastic scattering time. It is important that because of 
the destructive interferences between rotated skyrmions, 
the rate is suppressed to zero in a QSDpSS even at a finite 
temperature and the skyrmion charge Cu(i) is conserved 
as emphasised in a few occasions in the paper. 

In TSDpSSs, the Hopf term is irrelevant because L T <C 
A s ; following Eq. |Cl| , the partition function of a monopole 
configuration {x m } is 



lnG(l,l) cx L r , 



(70) 



Z = 



Nl 

N Q a =±l" a=l 



exp(- S%f + iXsQ«)- 
(68) 

And Sg 1 ~ A/T represents a thermal activation factor. 

Eq.[38] suggests that in 2d TSDpSSs, the skyrmions are 
condensed and < Q 2 >^ 0. The topological charge C;y(i) 
is not conserved. As emphasised before, when the parti- 
cle moves around a condensed skyrmion of an arbitrary 
size, it picks up a Berry's phase between zero and 2tt. 
The probability for a particle to propagate at a large dis- 
tance L is 



G(l,2)~£e*p(-tT m ),r = £; f dr-A(r,i); 



(69) 



A(i) is the vector potential due to the ith skyrmion and 
r m is the total phase acquired by the quasiparticle mov- 
ing along path m. Given that skyrmions are randomly 
nucleated in the 2D space and assuming the paths are 
cigar-like ones, we estimate 



with L m being the length of a characteristic loop m. The 
estimate shows the particle is localized because of the 
destructive interferences between different paths in the 
presence of thermally nucleated skyrmions. Or say it dif- 
ferently, a bare BdeG is infinitily massive because of the 
skyrmion cloud. 




Fig. 6 Quasiparticles in a cloud of skyrmions. The filled 
circles stand for the positive skyrmions and empty ones 
for negative skyrmions. The Berry's phase for a quasi- 
particle moving along a skyrmion could be between zero 
and 2n. 



For a similar reason, the expectation value in Eq,62 
which is equal to the Wilson loop integral of a gauge field, 
decays exponentially with an exponent proportional to 
PL T L2; P = exp(— Sq) is the thermal activation rate. 
The energy of a BdeG quasiparticle is proportional to the 
size of the system. And the bare BdeG particles could 
have been confined in a TSDpSS. In general, in a TS- 
DpSS, the mechanism of having zero energy skyrmions is 
absent and the Cym is not conserved. Bare BdeG quasi- 
particles interact with a linear potential and form bound 
states. 

Tab. 1 Comparison of quasiparticles 





Compos. 


Spin 


Qe.m. 


Qu(i) 


DoS 


Sta. 


SOpSS 


bare 
BdeG 


1/2 


0* 


i 


gapped 


Fermi 


QSDpSS 


BdeG+ 
Skyr. 


1/2 








gapped 


Bose 


TSDpSS 


BdeG+ 
BdeG 


0,1 








gapped 


Bosc 



*This is for the charge at e(k) = 0. Q e . m . , Qu(i) are 
charges defined with respect to fields A e . m . and A. 



C. Collective spin excitations 

As emphasised in the previous section, in TSDpSSs, 
the U(l) fields are confining. An excitation carrying a 
unit charge is confined and is forbidden. The collective 
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excitations have to be "charge neutral" with respect to 
A. A spin- 1/2 77 which carries a unit charge is absent 
in the excitation spectrum; only pairs of i] which carry 
zero U(l) charges exist as physical excitations in a form 
of the usual spin wave excitations. This is what happens 
in TSDpSSs. 

However, the topological term changes the nature of 
interactions mediated by A fields in QSDpSSs. In fact, 
an r\ quantum is bound with a skyrmion such that the 
bound state becomes a fermion. Each bosonic spin wave 
excitation is fractionalized into two elementary fermionic 
spinors hosted in skyrmions in the spin disordered limit. 
Each spinor-skyrmion composite carries a half spin but 
no U(l) charge, by contrast to a bare 77 excitation. 



(a) 




(b) 



Fig. 5 a) The confinement of 7r-disclinations(patterned 
circles) and half vortices (empty circles) in a SOpSS by a 
string, or a branch cut; b) the liberation of half-vortices 
in a SDpSS. 



Tab. 2 Comparison of collective spin excitations 





Compos. 


Spin 


Qe.m. 


Qu(i) 


DoS 


Sta. 


SOpSS 


Spin 

waves 


1 








gapless 


Bose 


QSDpSS 


Spinor + 
Skyr. 


1/2 








gapped* 


Fermi 


TSDpSS 


Spinor + 
Spinor 


0, 1 








gapped* 


Bose 



* However, it is unclear how the magnitude of the gap 
and the range of spin correlation are renormalized by 
the topological term. 



V. TOPOLOGICAL EXCITATIONS 



The SDpSS discussed here has the following order pa- 
rameters: 

< A >= 0, Tr < AA >^ 0, < exp(ix) >^ 0. (71) 

The existence of SC* with <A>=0,Tr<AA>^0, 
< exp(i%) >= 0, and other fractionalized states exam- 
ined in jll| is beyond the model studied here. Physically, 
the SDpSS has Josephson oscillations of 2eV frequency 
while in SC* the frequency is AeV. 

In the presence of spin-orbital couplings, the mean field 
solution indicates that the director of n points along ±e z 
direction and the internal space is [Z 2 x S 1 ]/Z 2 - How- 
ever, at the energy scale higher than the spin-orbit cou- 
pling ones, n would be free to rotate on a two-sphere. 
The spin order-disorder transition still takes place at a 
finite temperature below the superconductor-metal tran- 
sition temperature T c when the spin-orbit scattering rate 
is much smaller than Ao(0). As the spin is disordered, 
skyrmions are again zero energy objects and the discus- 
sions on the quasiparticles are still valid. 



VI. CONCLUSION 



For the sake of completeness, I am also listing some 
properties of vortices. The linear defects of a symme- 
try broken state with an internal space 1Z = [S 1 x S 2 ]/Z2 
have been recently discussed extensively in the context of 
Bose-Einstein condensates of 23 Na @5|]. In a SOpSS, 
the linear defects are superpositions of hc/Ae vortices and 
7r-disclinations because of the Z 2 symmetries in the prob- 
lem. An individual 7r-disclination or hc/Ae vortex carries 
a branch cut along which the phase jumps by 7r. The 
branch cut emitted by a hc/Ae vortex can end only at a tt- 
disclination, resulting in a linear potential between these 
objects and therefore confinement. And a bare hc/Ae 
vortex is forbidden because of the catastrophe of a cut. 
In a SDpSS, however, hc/Ae vortices can exist by its own 
right and are elementary excitations. 



In conclusion, we also would like to remark that some 
aspects of the quasiparticles in spin disordered supercon- 
ductors considered here reminisce the chiral-bag defect 
model for the nucleoli |34|6q]. The presence of spin- 1/2 
bosonic chargeless BdeG excitations in a SDpSS is an 
example of fermi number fractionalization; it belongs to 
the same class phenomenon as the mid-gap quasiparti- 
cles hosted in domain wall excitations in one dimension 
polyacetylene || . and the statistical transmutations pro- 
posed in some magnetic models [ pOp^p^ ]. Though in 
this case all the excitations in the spin rotation invariant 
superconducting state remain gapped, the spectral prop- 
erty should be very different from that in conventional 
spin triplet superconductors when the S 2 symmetry is 
broken |66| . 
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APPENDIX A: DERIVATION OF EQ.|l3j USING 
SUB-LATTICES 

Eq.|l3] can also be obtained by considering two sublat- 
tices A and B with u(JTa) and t>(fii 2 ) labeling the state 
at site i of the lattice A and B. At each site the total 
spin of A and B is one. One can use the results for u(Q) 
in Eq.5 to obtain the effective field 



F„ 



F A 



F B 



p A,B 



— 2 ta Pl Si l,2 



dn'l 2 dW l2 



(Al) 



Taking into account f2i j2 = n — ^, one obtains the re- 
sults previously derived in Eq.|l^. From here it is clear 
how the contributions from u(fli) and u(J~2 2 ) cancel with 
each other. 



APPENDIX B: DERIVATION OF 
COEFFICIENTS IN EQ 



+ 00 N , 

z = e ™ n c e 



iV=0, Q = l 



Qc=±l 



exp (i^2jB{Q a ,T a ) + i^2xsQa - Sq), 



QaQfi 



^ + (CD 



2 5 (|r^-r Q |2 + (^-r Q ) 2 ) V2 ' ^ 2 9 



Note that the first phase factor in the exponent arises 
from the Hopf term, distinguishing this system from the 
usual compact U{1) theory. Xs is introduced as a Sine- 
Gordon field coupled to a monopole charge Q a . 



APPENDIX D: GAUGE FIELDS DURING THE 
GAUGE TRANSFORMATION IN EQ.|| 

The gauge fields in between (0 < a < 1) can be choosen 
to be 



T-Tl 



T 2 ~T 



2(r 2 + (r _ Tl) 2 )3 /2 2(r2 + (r _ T2)2) 3/2 

H x = [x cos (3{t) +y sin 0(t)] x 

1 1 



'2(r» + (r-7i)»)8/» 2 ( r 2 + (T _ T2) 2)3/2 
Hy = [— x sin (3(t) + ycos/3(r)] x 



1 



1 



L 2(r2 + ( Tl -T)2)3/2 2(r2 + (T „ T2)2) 



3/2 J 



(Dl) 



7rfcT ^ AQV a v 6 



r 



E 



£?(k)s 



A 2 
^0 



V 



k ^ 2 +J E(k) 2 } 



(Bl) 



and /3(ti) — /3(r 2 ) — n2n. T\(a) < T 2 (a) varies from 
0(a = 0) to TL T /2 (a = 1/2) and back to at a = 1. In 
this configuration, the skyrmion charge varies as a func- 
tion of r, 



C U (i) = l[Q(r-T 1 )-e(r 1 -r) 
+9(t 2 - r) - 6(r - r 2 )]. 



(D2) 



(2n + lWfcT, n = 0, 1,2... are the Masubara fre- 



quencies. Eq.Bl leads to Eq.|23[ 
iq 2 in Eq.|23| can be calculated as 

d 2Aq(T) jg(k) 
^i(k,r)-- ae2(k) f;(k) tan(— ) 



F 2 (k,T) = - 



A 2 



5 



£ 3 (k) de(k) v £(k) 



+ — ^ 



(B2) 



APPENDIX C: THE ACTION FOR SPACE-TIME 
MONOPOLES 

The action for space-time monopoles is 
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